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Abstract 

There has been a great deal of interest in understanding the properties of quantum chromo- 
dynamics (QCD) for a finite value of the chemical potential and for finite temperature. Studies 
have been made of the restoration of chiral symmetry in matter and at finite temperature. The 
phenomenon of deconfinement is also of great interest, with studies of the temperature depen- 
dence of the confining interaction reported recently. In the present work we study the change of 
the properties of light mesons as the temperature is increased. For this study we make use of 
a Nambu-Jona-Lasinio (NJL) model that has been generalized to include a covariant model of 
confinement. The parameters of the confining interaction are made temperature-dependent to take 
into account what has been learned in lattice simulations of QCD at finite temperature. The con- 
stituent quark masses are calculated at finite temperature using the NJL model. A novel feature 
of our work is the introduction of a temperature dependence of the NJL interaction parameters. 
(This is a purely phenomenological feature of our model, which we do not attempt to derive from 
more fundamental considerations.) With the three temperature-dependent aspects of the model 
mentioned above, we find that the mesons we study are no longer bound when the temperature 
reaches the critical temperature, Tc, which we take to be 170 MeV. We believe that ours is the first 
model that is able to describe the interplay of chiral symmetry restoration and deconfinement for 
mesons at finite temperature. The introduction of temperature-dependent coupling constants is a 
feature of our work whose further consequences should be explored in future work. 
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I. INTRODUCTION 



In recent years we have studied a generalized version of the Nambu-Jona-Lasinio (NJL) 
model which includes a covariant model of Lorentz- vector confinement [1-5]. Extensive 
applications have been made in the study of light mesons, with particularly satisfactory 
results for the properties of the r/(547), rj' {958) and their radial excitations [6]. Since the 
modifications of the confining potential at finite temperature have recently been obtained 
in lattice simulations of QCD [7-12] (see Fig. 1), we became interested in introducing that 
feature in our generalized NJL model, whose Lagrangian is 

i=0 

+^{det[g(l + 75)g] + det[g(l - 75)?]} 

+^conf ■ (1-1) 

Here, the A*(i = 0, • • • , 8) are the Gell-Mann matrices, with A° — -\/2/3 1, m° = 
diag (m°, m^, m°) is a matrix of current quark masses and Cconf denotes our model of con- 
finement. 

The dependence of the constituent masses of the NJL model upon the density of nuclear 
matter has been studied in earlier work, where we also studied the deconfinement of light 
mesons with increasing matter density [13]. In that work we introduced some density depen- 
dence of the coupling constants of the model. However, that was not done in a systematic 
fashion. As we will see, in our study of the temperature dependence of the constituent 
masses, as well as in the calculation of temperature-dependent hadron masses, we introduce 
a temperature dependence of the coupling constants, which for T — Tc, represents a 17% 
reduction of the magnitude of these constants. Since the use of temperature-dependent 
coupling constants is a novel feature of our work, we provide some evidence that such tem- 
perature dependence is necessary to create a formalism that is consistent with what is known 
concerning QCD thermodynamics. This aspect of our work is discussed in the Appendix. 

The organization of our work is as follows. In Section II we review our treatment of 
confinement in our generalized NJL model and describe the modification we introduce to 
specify the temperature dependence of our confining interaction. In Section III we discuss 
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FIG. 1: A comparison of quenched (open symbols) and unquenched (filled symbols) results for the 
interquark potential at finite temperature [7]. The dotted line is the zero temperature quenched 
potential. Here, the symbols for T = O.SOTc [open triangle], T = 0.88Tc [open circle], T = 0.94Tc 
[open square], represent the quenched results. The results with dynamical fermions are given 
at T = 0.68Tc [solid downward-pointing triangle], T = O.SOTc [solid upward-pointing triangle], 
T = O.BSTc [solid circle], and T = OMTc [solid square]. 

the calculation of the temperature dependence of the constituent masses of the up, down 
and strange quarks. In Section IV we describe how random-phase-approximation (RPA) 
calculations may be made to obtain the wave function amplitudes and masses of the light 
mesons. (We consider the properties of the tt, K, Jq, Cq and Kq mesons and their radial 
excitations in this work.) In Section V we provide the results of our calculations of the 
meson spectra at finite temperature. Finally, Section VI contains some further discussion 
and conclusions. 
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II. COVARIANT LORENTZ- VECTOR MODEL OF CONFINEMENT 



We have published a number of papers in which we have described our model of 
confinement [1-6]. In this work we provide a review of the important features of the model. 
As a first step, we introduce a vertex function for the confining interaction. [See Fig. 2.] 
For example, the vertex useful in the study of pseudoscalar mesons satisfies an equation of 
the form [6] 

dV 



[rSa{P/2 + k') 



xr,,ab{P, k')S,{-P/2 + k')^^]V^{k - k') , 



(2.1) 



where Sa{P/2 + k') = [f/2 + I/t' — rria + irj] ^ . Here, and are the constituent quark 
masses. We consider the form V'"{r) = «;rexp[— /xr] and obtain the Fourier transform 



V^{k-k') = -Sttk 



1/' 



(2.2) 



[{k-k'f + n^]^ [{k-k')^ + H^f 

which is used in Eq. (2.1). Note that the matrix form of the confining interaction is 
V'-'{k — k') = 7^(l)y'"(/c — A;')7^(2), since we are using Lorentz-vector confinement. The 
form of the potential may be made covariant by introducing the four-vectors 

{k ■ P)P'' 



k^^k^- 



p2 



and 



k'l" = k'l" - 



{k' ■ P)Pi' 

P2 



(2.3) 



(2.4) 



We then define 



4/^2 



[-{k-k')^ + 11^]^ [-{k-k')^ + li^]^ 

which reduces to V'"\k — k') of Eq. (2.2) in the meson rest frame, where P = 0. 

It is also useful to define scalar functions F^ ~(,(P, k) and F5^^(P, k) [6]. We introduce 



(2.5) 



2m„ 



and 



" ^ ^ 2m„ ' 



(2.6) 



(2.7) 



P/2 + k P/2 + k p/2 + k' P/2 + k 

->• r >- 





-P/2 + k -P/2 + k -P/2 + k' -P/2 + k 



FIG. 2: The figure provides a schematic representation of Eq. (2.1) for the confining vertex function 

where A;'^ = [Ea{k), k] and /c'^ = [—Ea{k),k], and define 

M/\mAP^ k)Ai-\-k) = r+- (P, k)A<^^\kh,Ai-\-k) , (2.8) 

and 

Ai-\-mAP, k)A^^\k) = r,;jP, k)Ai-\-kh,A^;\k) . (2.9) 

We have obtained equations for r^~(P, k) and r^"^(P, k). For example, with = m^, we 
have 



rr{P',\k\)^i-J- 



-2E{k)E{k') 



Tt-iP'Ak'DV^ik-k') 

P^-2E{k') 



E{k)E{k') 

Note that r^~(P°, |^'|) = 0, when P° = 2E(k'), so that one need not introduce a small 
quantity, ie, in the denominator on the right-hand side of Eq. (2.10). (Alternatively, we 
may note that V^^i^P^, \ k\) = when the quark and antiquark in Fig. 2 are on mass shell.) 
The confinement vertex functions defined in our work may be used to calculate vacuum 
polarization functions which are real functions. The unitarity cut, that would otherwise 
be present, is eliminated by the vertex functions which vanish when both the quark and 
antiquark go on mass shell [1-6]. 

We have presented Eq. (2.10), since we wish to stress that for the study of hght mesons 
the constituent quark mass, m, is of the order of \k\, so that \k\/m is not small. Therefore, 
the term in the square bracket on the right-hand side of Eq. (2.10), which would be equal to 
—1 in the nonrelativistic limit, provides quite important relativistic corrections to the form 
given in Eq. (2.5). 

In Fig. 3a we show the homogeneous equation for the confining vertex. The solution of 
the homogeneous equation allows us to construct the wave functions bound in the confining 



p/2 + k p/2 + k' 




P/2 + k 

~> — ^ — > ^ — > — ■ 




-< — <- 



-P/2+k , , -P/2 + k' -P/2 + k 

(a) 



P/2 + k P/2 + k' P/2 + k' P/2 + k* P/2 + k 

-P/2 + k -P/2 + k' -P/2 + k' -P/2 + k' -P/2 + k 

(b) 



FIG. 3: a) Bound states in the confining field (wavy fine) may be found by solving the equation 
for the vertex shown in this figure, b) Effects of both the confining field and the short-range NJL 
interaction (filled circle) are included when solving for the vertex shown in this figure. 



field. For example, we may define 

- 2E{k) 

and 



= '^f'l'S (2.11) 



which we may term the "large" and "small" components of the wave function of the bound 
state with energy Pf. In Fig. 3b we show the equation for the vertex function that includes 
both the effects of the short-range NJL interaction and the confining interaction. We will 
return to a consideration of the equations obtained in an analysis of Fig. 3b in Section IV 
where we consider the RPA equations of our model. 

In order to motivate our treatment of the temperature dependence of the confining in- 
teraction, we have presented some results obtained with dynamical quarks (filled symbols) 
in Fig. 1. The fact that the potential becomes (approximately) constant for r > 1 fm is 
ascribed to "string breaking" in the presence of dynamical quarks. (Note that, upon string 
breaking, the force between the infinitely massive quark and antiquark vanishes.) 
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For our calculations we have used // = 0.010 GeV and k — 0.055 GeV^ in the past. In 
order to introduce temperature dependence in our model, we replace V'^{r) — «;rexp[— //r] 
by V^'"(r, T) — «;r exp[— //(T)r], with 

"t^)- 1-0.7" W - 

where jiQ — 0.01 GeV. Values of V'-^{r,T) for various values of the ratio T/Tc are given in 
Fig. 4. We remark that, while V'-^{r) — > for large r, the bound-state solutions found for 
l/'-^(r) are largely unaffected, since barrier penetration effects are extremely small in our 
model. The maximum value of the potential is V^^^^ = K/fie with the corresponding value 
of Tjnax = V/"- Thus, in the study of the bound states, our model is essentially equivalent 
to one with V'~^{r) — Kr exp[—fj,r] for r < r^ax a-nd \^*^(r) = Vmax ^ ''''max- The same 
remarks pertain, if we replace ji by //(T") of Eq. (2.13). With that replacement, we have 

Wr)^ "'^-"-^™ . (2.14) 

We note that Vmax is finite at T = Tc , a result that is in general accord with what is found 
in lattice calculations of the interquark potential for massive quarks. 



III. CALCULATION OF CONSTITUENT QUARK MASSES AT FINITE TEM- 
PERATURE 

In an earlier work we carried out a Euclidean-space calculation of the up, down, and 
strange constituent quark masses taking into account the 't Hooft interaction and our con- 
fining interaction [14]. The 't Hooft interaction plays only a minor role in the coupling of 
the equations for the constituent masses. If we neglect the confining interaction and the 't 
Hooft interaction in the mean field calculation of the constituent masses, we can compensate 
for their absence by making a modest change in the value of Gg- [See Eq. (1.1).] For the 
calculations of this work we calculate the meson masses using the formalism presented in 
the Klevansky review [15]. (Note that our value of Gs is twice the value of G used in that 
review.) The relevant equation is Eq. (5.38) of Ref. [15]. Here, we put // = and write 

m{T) = m° + AGNc"^ [ dp^ tanh(^/3£:p) , (3.1) 

TT^ Jq hp 1 
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FIG. 4: The potential V'^{r,T) is shown for T/Tc = [sohd line], T/Tc = 0.4 [dotted line], 
T/Tc = 0.6 [dashed line], T/Tc = 0.8 [dash-dot line], T/Tc = 0.9 [short dashes], T/Tc = 1-0 
[dash-dot-dot line]. Here, F^(r,T) = /cr exp[-/x(r)r], with /x(r) = 0.01GeV/[l - 0.7(r/rc)2]. 

where A = 0.631 GeV is a cutoff for the momentum integral, (3 — 1/T and Ep — [p^ + 
m^(T)]^/^. In our calculations we replace G by Gs{T)/2 and solve the equation 



with Gs{T) = 11.38[1 - 0.17 T/TJ GeV, m° = 0.0055 GeV and m° = 0.120 GeV. Thus, we 
see that Gs{T) is reduced from the value Gs{0) by 17% when T — Tc. The results obtained 
in this manner for mu{T) and ms{T) are shown in Fig. 5. Here, the temperature dependence 
we have introduced for Gs(T) serves to provide a somewhat more rapid restoration of chiral 

symmetry than that which is found for a constant value of Gs- That feature and the 
temperature dependence of the confining field leads to the deconfinement of the light mesons 
considered here at T = Tc. (See Section V.) 



m{T) 



m° + 2Gs{T)Nc 



m 




(3.2) 
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FIG. 5: Temperature dependent constituent mass values, mu{T) and ms{T), calculated using Eq. 
(3.2) are shown. Here m° = 0.0055 GeV, m° = 0.120 GeV, and G{T) = 5.691 [1 - 0.17(r/rc)], if 
we use Klevansky's notation [15]. (The value of Gs of Eq. (1.1) is twice the value of G used in 
Ref. [15]). 

IV. RANDOM PHASE APPROXIMATION CALCULATIONS FOR MESON 
MASSES AT FINITE TEMPERATURE 

The analysis of the diagrams of Fig. 3b gives rise to a set of equations for various 
vertex functions. These equations are of the form of relativistic random-phase-approximation 
equations. The derivation of these equations for pseudoscalar mesons is given in Ref. [6], 
where we discuss the equations for pionic, kaonic and eta mesons. The equations for the 
eta mesons are the most comphcated, since we consider singlet-octet mixing as well as 
pseudoscalar-axial-vector mixing. In that case there are eight vertex functions to consider, 
FpQ, Tpg, r^g, TpQ, r^Q, Fpg, r^g, where P refers to the 75 vertex and A refers to 
the 7o75 vertex, which mixes with the 75 vertex. Corresponding to the eight vertex functions 
one may define eight wave function amplitudes [6]. 

The simplest example of our RPA calculations is that of the oq mesons, where there are 
only two vertex functions F^~ and F^"*" to be calculated [4]. Associated with F^~ and F^"*" 
are two wave functions (pg and 0^, which are the large and small components respectively. 
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In vacuum one has coupled equations for these wave functions. 

2E^{k)(l)+{k) + J dk' [Hc{k,k') + HNjL{k,k')](l)+{k') (4.1) 
+ Jdk' H^jL{k,k')r{k') = P'<P+{k) , 

- 2E^{k)<p-{k) - j dk'[Hc{k, k') + HNjL{k, k')]<p-{k') (4.2) 
" dk' HnjL{Kk')(t)-^{k') ^ P''(t)-{k) , 



where E^(k) = f^ + m^Ji/s^ 



1 [2V,^{k,k')k'k'' + mlkk'Vf{k,k')] 

Hcik,k ) = £;„(fc)E„(fcO ' ^'-'^ 



and 



^"-^^ - (2^0^ ^SSM^O ■ ^ 

In Eq. (4.4), is the effective couphng constant for the Qq mesons, which depends upon 
the values of Gs, Gd and the vacuum condensates. These relations for the various coupling 
constants, Ga^i Gtt, Gki Gf^, Gk*, etc. may be found in Ref. [16]. In Eq. (4.3) we have 
introduced 

Vf^{k,k') = J dxPi{x)V^{k-k'). (4.5) 

Here, x = cos9 and Pi{x) is a Legendre function. The terms exp[— /c^/2q;^] and exp[— /c'^/2q;^] 
are regulators with a = 0.605 GeV. 

In order to solve these equations at finite temperature, we replace m„, Gao and /iq by 
mu(T), Gao{T) and //(T"). The values of m{T) are given in Fig. 5, and we recall that 
H{T) = no/[l - 0.7(r/rc)2]. In the RPA, the solutions of Eqs. (4.1) and (4.2) come in 
pairs. For a state of energy P° there is another state with energy —Pi- Since the RPA 
Hamiltonian is not Hermitian, it is possible to obtain imaginary values for the energy. That 
is a signal of the instability of the ground state of the theory and requires that the problem be 
reformulated to obtain a stable ground state. This problem does not arise in the calculations 
reported in this work. In particular, the use of the temperature-dependent values of Gt^{T) 
avoids the appearance of pion condensation in the formalism. 
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V. RESULTS OF NUMERICAL CALCULATIONS OF MESON MASSES AT 
FINITE TEMPERATURE 

As noted earlier, the RPA equations for the Oq mesons are given in Ref. [4] and those 
for the TT, K and rj mesons are given in Ref. [6]. The equations needed in the study of 
the /o mesons are to be found in Ref. [5], while the RPA equations for the study of the 
Kq mesons are to be found in the Appendix of Ref. [3]. In Figs. 6-10 wc present our 
results for the tt, K, qq, /q and Kq mesons. The values of the coupling constants used are 
given in the figure captions. The reduction of the number of bound states with increasing 
temperature can be understood by noting that for the tt and Oq mesons the continuum of 
the model hes above Vmax{T) + 2m„(T), while for the K and Kg mesons the continuum 
lies above VmaxiT) + mu(T) + ms{T). The situation is more complex in the case of the /o 
mesons which contain both ss, uu and dd components. The bound states lie below the ss 
continuum which begins at VmaxiT) + 2ms{T). However, we note that the absence of bound 
states a,t T — Tc for all the mesons considered here is due to the reduction of the value of 
the confining potential and of the constituent quark masses. 

It is of interest to note that the mass values of the Oq and /o mesons tend toward de- 
generacy with the pion as T ^ Tc . However, the mesons disassociate before a greater 
degeneracy is achieved. That is in contrast to the results obtained in the SU(2) formalism 
considered by Hatsuda and Kunihiro [16]. Since these authors do not include a model of 
the confinement-deconfinement transition, they are able to see the approximate degeneracy 
of the sigma meson and the pion with increasing temperature. It is also worth noting that 
in the SU(3) formahsm the sigma meson is replaced by the /o(980) as the chiral partner of 
the pion. 

In order to demonstrate the interplay of chiral symmetry restoration and dissolution of 
our meson states, we have performed calculations in which the quark masses are unchanged 
from their value at T = 0. For the pion, with ruu — ma — 0.364 GeV and for T/Tc — 0.95, 
we find bound states at 0.530, 1.242 and 1.305 GeV. If we also consider the values of the 
coupling constants and quark masses fixed at their T — values, bound states are found 
at 0.102, 2.248 and 1.298 GeV when T/T^ = 0.95. A similar analysis for the kaon states 
yields bound states at T = Tc of 0.738, 1.395 and 1.444 GeV, if we put m„ = 0.364 GeV and 
rUs = 0.565 GeV. If, in addition, we neglect the temperature dependence of the coupling 
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FIG. 6: The mass values of the pionic states calculated in this work with Gt^{T) = 13.49[1 — 
O.ITT/Tc] GeV, Gv{T) = 11.46[1 - O.ITT/Tc] GeV, and the quark mass values given in Fig. 
5. The value of the pion mass is 0.223 GeV at T/Tc = 0.90, where m„(T) = 0.102 GeV and 
mgiT) = 0.449 GeV. The pion is bound up to T/Tc = 0.94, but is absent beyond that value. 

constants, we have bound kaon states at 0.482, 1.440 and 1.439 GeV. 

In the case of the states, putting m„ = = 0.364 GeV yields a single state at 0.960 
GeV at T = Tc , if we neglect the temperature dependence of the coupling constant. If we 
maintain the temperature dependence of the coupling constant, we find a single state at 
1.067 GeV, if the quark masses are kept at their T = values. 

Prom this analysis, we see that the reduction of the quark masses with increasing tem- 
perature, which represents a partial restoration of chiral symmetry, is an essential feature of 
our model. In the past, lattice simulations of QCD have indicated that deconfinement and 
restoration of chiral symmetry take place at the same temperature. Since our Lagrangian 
contains current quark masses for the up, down and strange quarks, chiral symmetry is not 
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FIG. 7: Mass values of kaonic states calculated with Gk{T) = 13.07[1 - O.ITT/Tc] GeV, Gv{T) = 
11.46[1 — 0.17 r/Tc] GeV, and the quark mass values given in Fig. 5. The value of the kaon mass 
is 0.598 GeV at T/Tc = 0.95, where m„(r) = 0.075 GeV and ms{T) = 0.439 GeV. 

completely restored at the higher temperatures in our model. However, the model does 
exhibit quite significant reductions of the up and down constituent quark masses for T > Tc 
(See Fig. 5), so that deconfinement is here associated with a significant restoration of chiral 
symmetry. 



VI. DISCUSSION 



In recent years we have seen extensive applications of the NJL model in the study of 
matter at high density [17-21]. There is great interest in the diquark condensates and color 
superconductivity predicted by the NJL model and closely related models that are based 
upon instanton dynamics. It is noted by workers in this field that the NJL model does not 
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FIG. 8: Mass values for the oq mesons calculated with Gao{T) = 13.1[1 - 0.17 T/Td GeV, and the 
quark mass values given in Fig. 5. The value of the oq mass at T/Tc = 0.95 is 0.416 GeV. 

describe confinement, with the consequence that one can not present a proper description 
of the hadronic phase that exists at the smaller values of the temperature and density in 
the QCD phase diagram. Thus, the attitude adopted is that, if one works in the deconfined 
phase, the NJL model may provide a satisfactory description of the quark interaction. In the 
present work we have modified the NJL model so that we can describe light mesons and their 
radial excitations, as well as the confinement-deconfinement transition at T = Tc . In an 
earlier work we studied the confinement-deconfinement transition for finite matter density at 
T = [13]. (A more comprehensive study would include both finite temperature and finite 
density.) As in the present work, in which we introduced a temperature-dependent coupling 
constants, we used density- dependent coupling constants in Ref. [13]. If such dependence 
exists, it would have important consequences for the study of dense matter using the NJL 
model. 




14 
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FIG. 9: Mass values of the /o mesons calculated with Goo{T) = U.25[l -0. 17 T/Tc] GeV, GssiT) = 
10.65[1 - O.ITT/Tc] GeV, GosiT) = 0.495[1 - O-lTT/rJ GeV, and Gso{T) = GosiT) in a singlet- 
octet representation. The quark mass values used are shown in Fig. 5. The fo has a mass of 0.400 
GeV at T/Tc = 0.95. 

One interesting feature of our results is that both the lowest Gq and fo states move toward 
degeneracy with the pion as the temperature is increased. However, the system is deconfined 
before such degeneracy can be exhibited. 

We stress that the restoration of chiral symmetry is intimately connected with the dis- 
solution of our meson states at T = Tc . As we saw in the discussion toward the end of 
Section V, the various mesons studied here still have bound states at T = Tc , if only the 
temperature dependence of the confining field is included in the model. 

The behavior of charmonium across the deconfinement transition has recently been stud- 
ied using lattice simulations of QCD [22]. The authors of Ref. [22] point out that, unlike 
the case of light mesons, charmonia may exist as bound states even after the deconfinement 
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FIG. 10: Mass values obtained for the K* mesons with Gk*{T) = 10.25[1 - 0.17 T/Tc] GeV, and 
the quark mass values shown in Fig. 5. 

transition. They state: "Our studies support the sequential pattern for charmonium disso- 
lution obtained from potential model studies, where the broader bound states (the scalar 
and axial vector channels) dissolve before the pseudoscalar and vector channels [23]. The 
pseudoscalar and vector channels are seen to survive as bound states still at 1.25Tc and 
probably dissolve after 1.5Tc." 

APPENDIX 

Since the introduction of temperature-dependent coupling constants for the NJL model 
is a novel feature of our work, we provide arguments in this Appendix to justify their 
introduction. We make reference to Fig. 1.3 of Ref. [24]. That figure shows the behavior of 
the ratio e/T^ and 3P/T^ for the pure gauge sector of QCD. Here e is the energy density 
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and P is the pressure. Ideal gas behavior imphes e — ?>P. The values of e/T^ and 3P/T'^ 
are compared to the value ese/T^ = 87r^/15 for an ideal gluon gas. It may be seen from the 
figure that at T = ST^ there are still significant differences from the ideal gluon gas result. 
Deviations from ideal gas behavior become progressively smaller with increasing T/T^ and 
could be considered to be relatively unimportant for T/T^ > 5. 

To provide evidence for temperature-dependent coupling constants, we discuss the cal- 
culation of hadronic current correlators in the deconfined phase. The procedure we adopt 
is based upon the real-time finite-temperature formalism, in which the imaginary part of 
the polarization function may be calculated. Then, the real part of the function is obtained 
using a dispersion relation. The result we need for this work has been already given in the 
work of Kobes and SemenofF [25]. (In Ref. [25] the quark momentum is k and the antiquark 
momentum is k — P. We will adopt that notation in this section for ease of reference to 
the results presented in Ref. [25].) With reference to Eq. (5.4) of Ref. [25], we write the 
imaginary part of the scalar polarization function as 

{(1 - ni(k) - n2(A;))5(p° - E,(k) - E^ik)) 
-Mk) - n^mSip'' + E^{k) - E^m 
-{n2{k) - ni(A;))5(p° - Ei{k) + E^ik)) 
-(1 - nr{k) - n2ik))5{p' + E,{k) + E^ik))} . 

Here, Ei{k) = [k'^ + m^(r)]V2. Relative to Eq. (5.4) of Ref. [25], we have changed 
the sign, removed a factor of and have included a statistical factor of 2Nc, where the 
factor of 2 arises from the flavor trace. In addition, we have included a Gaussian regulator, 
exp[— ^^/a^], with a — 0.605 GeV, which is the same as that used in most of our applications 
of the NJL model in the calculation of meson properties. We also note that 

Mk) - ^pj,) ^ ^ , (A2) 

and 

^2(^) = e/^^^(fe) + r ^^^^ 

For the calculation of the imaginary part of the polarization function, we may put k^ — 
mf (T) and {k—pY = m|(T), since in that calculation the quark and antiquark are on-mass- 
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shell. We will first remark upon the calculation of scalar correlators [25] . In that case, the 
factor 13s in Eq. (Al) arises from a trace involving Dirac matrices, such that 

13s = -Tr[(^ + ^,)(^ _ + m2)] (A4) 
= 2p2-2(ml + m2)^ (A5) 

where mi and m2 depend upon temperature. In the frame where P = 0, and in the case 
mi — 1712, we have Ps — 2P|(1 — A^m^ /Pq). For the scalar case, with mi —m2, we find 

Im J-(P^T) - ^ (l - - , (A6) 

where 

= ^-m\T). (A7) 

We may evaluate Eq. (2.8) for m{T) = mu{T) = md{T) and define Im J^(P^,T). Then 
we put m(T) = ■ms{T), we define Im Jf(P^,T). These two functions are needed for a 
calculation of the scalar-isoscalar correlator. The real parts of the functions J^(P^,T) and 
Jf (P^,T) may be obtained using a dispersion relation, as noted earlier. 

For pseudoscalar mesons, we replace Ps by 

/3p = -TV[i75(|^ + mi)i75(|^ - f + ms)] (A8) 
= 2P2 - 2(mi - ms)^ , (A9) 

which for mi = m2 is Pp — 2Pq in the frame where P — 0. We find, for the tt mesons, 

Im J''(P^T) = ^ (l - '-^y\-'^''"ll - 2u,m . (AlO) 

where k^ = Pq/A — m'^{T), as above. Thus, we see that, relative to the scalar case, the 
phase space factor has an exponent of 1/2 corresponding to a s-wave amplitude, rather than 
the p-wave amplitude of scalar mesons. For the scalars, the exponent of the phase-space 
factor is 3/2, as seen in Eq. (A6). 

For a study of vector mesons we consider 

P^. = T^[7m(^ + -r + 1712)] , (All) 

and calculate 

gf^^P^, = 4[P2 -ml-ml + Am,m2] , (A12) 
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which, in the equal- mass case, is equal to 4Pq^ + 8m^(T), when mi = m2 and P — 0. Note 
that for the elevated temperatures considered in this work m„(T) = rad{T) is quite small, 
so that APq + 8m^(T) can be approximated by AP^ when we consider the p meson. 

We now consider the calculation of temperature-dependent hadronic current correlation 
functions. The general form of the correlator is a transform of a time-ordered product of 
currents, 

C(P^T) = i j d^xe'P-'' « T{j{x)j{0) » , (A13) 

where the double bracket is a reminder that wc arc considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp,i{x) = 
q{x)ij5X'q{x), where, in the case of the tt mesons, i — 1,2, and 3. For the study of 
pseudoscalar-isoscalar mesons, we again introduce jp,i{x) — q{x)X^q{x), but here i — 
for the flavor-singlet current and i = 8 for the flavor-octet current. 

In the case of the tt mesons, the correlator may be expressed in terms of the basic vacuum 
polarization function of the NJL model, Jp{P^,T). Thus, 

where Gt^{T) is the coupling constant appropriate for our study of the tt mesons. We have 
found G7r(0) = 13.49 GeV~^ by fitting the pion mass in a calculation made at T = 0, with 
= i^d = 0.364 GeV. 

The calculation of the correlator for pseudoscalar-isoscalar states is more complex, since 
there are both flavor-singlet and flavor-octet states to consider. We may define polarization 
functions for u, d and s quarks: Ju{P'^,T), Ja{P'^,T) and Js{P'^,T). In terms of these 
polarization functions we may then define 

Joo{P\T) = ^[J„(P^^) + Jd{P\T) + Js{P\T)] , (A15) 

MP',T) = ^[j„(p^r) + Jd(p^^) - 2Js{p^t)] , (Aie) 



and 



J88(p^T) = ^[J„(p^T) + j,(p2,T) +4J,(p^T)] . (ai?) 
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We also introduce the matrices 



J{P'^,T) 



MP\T) Jo8{P',T) 



(A18) 



G{T) 



Goo{T) Gos{T) 
Gso{T) Gss{T) 



(A19) 



and 



C{P^,T) 



Coo(P^T) Co8(P^T) 

C8o(P^T) CssiP^T) 



(A20) 



We then write the matrix relation 



G{P^,T) = J{P^,T)[1-G{T)J{P^,T)] 



-1 



(A21) 



The use of our energy-dependent coupling constants is meant to be consistent with the 

approach to asymptotic freedom at high temperature. In order to understand this feature in 
our model, we can calculate the correlator with constant values of Gqq, Ggs S'-iid Cqs and also 



with Goo{T) = Goo[l - O.nT/T^], etc. (In this work we use Cqo = 8.09 GeV-^, Ggs = 13.02 
GeV-2 and Gqs = -0.4953 GeV-^) 



We now consider the values of ImC88(-P^) for T/Tc — 4.0. In Fig. 11 we show the values 
of ImC88(i^'^) calculated in our model with temperature-dependent coupling parameters as 
a dashed line. The dotted line shows the values of the correlator for Goo = Gss = Gqs = 0, 
while the solid line shows the values when the coupling parameters are kept at their values 
at T = 0. We see that we have resonant behavior in the case the parameters are temperature 
independent. 

In Fig. 12 we show similar results for T/Tc — 5.88. Here the temperature-dependent 
coupling constants are equal to zero, so that the lines corresponding to the dashed and 
dotted lines of Fig. 11 coincide. The solid line again shows resonant behavior at a value 
of T/Tc where we expect only weak interactions associated with asymptotic freedom. We 
conclude that the model with constant values of the coupling parameters yields unacceptable 
results, while our model, which has temperature-dependent coupling parameters, behaves as 
one may expect, when the results of lattice simulations of QCD thermodynamics are taken 
into account. 
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FIG. 11: The imaginary part of the correlator CssiP'^) is shown for T/Tc = 4.0. The dashed hne 

is the result for the temperature-dependent coupling parameters of our model, while the solid line 

represents the results for coupling parameters kept at their T = values. The dotted line shows 

the values of the correlator when the coupling parameters are set equal to zero. 
0.10 n 




P' (GeV') 



FIG. 12: The imaginary part of the correlator Css^P'^) is shown for T/Tc = 5.88. [See caption to 
Fig. 11.] Here the dashed and dotted lines of Fig. 11 coincide. 
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